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A method of processing measurement while drilling data. 

@ A method and system are disclosed for the 
transformation of measurement while drilling 
(MWD) data (Fi) acquired during various time 
intervals into corresponding equal depth inter- 
vals. Once the data are transformed into the 
depth domain, depth based statistical, filtering 
sensor resolution matching and depth shifting 
techniques are utilized. The transformation pro- 
cess maximizes sensor vertical resolution while 
minimizing observed and statistical errors as- 
sociated with sensor response. 
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The present invention relates to processing of data obtained during the drilling of a well borehole. These 
data are typically produced by sensors positioned in the vicinity of the drill bit and are telemetered to the surface 
In real time for processing, or alternatively recorded down hole in real time for subsequent processing and ana- 
lysis. 

5 In the exploration and production of hydrocarbons, numerous techniques are used to analyze earth for- 

mations penetrated by a borehole. Wireline measurements involve a sensor bearing logging tool or sonde 
which is lowered and raised at a constant predetermined rate within the borehole. Data from the downhole sen- 
sors are transmitted to the surface by way of the cable used to raise and lower the logging tool. The resulting 
"logs" are graphical or digital tabulations of the measured parameters of interest as a function of depth within 

10 the well bore and are used to analyze prospective formations for hydrocarbon content. Data in wireline logging 
operations are typically collected over evenly spaced depth intervals since the final product is a depiction of 
measured parameters as a function of depth. Collection of data in this manner is controlled by depth indicator 
signals generated by the system as the speed at which the logging tool Is raised or lowered within the well 
bore. Wireline measurements must be made after the well bore has been drilled. Such measurements can be 

15 adversely affected by the build up of mud cake from the drilling fluid and/or other physical changes which occur 
in the vicinity of the borehole in the time interval between the drilling and logging operation. 

There are many advantages in making parametric measurements of interest during the actual drilling op- 
eratton. The formation is relatively unperturbed and valuable information concerning the petrophysics of the 
penetrated formations and the physk:al condition of the borehoie can be obtained in real time or near real time. 

20 Sensors can indicate in real time when the geological formation of interest has been reached. High pressured 
formations can be detected at the time of penetratk>n thus permitting steps to be taken to prevent damage to 
the well bore, drilling equipment and personnel. Measurement while drilling (MWD) is accomplished by placing 
sensors in the vicinity of the drill bit and either (a) telemetering data to the surface in real time or (b) recording 
the data downhole to be retrieved and processed intermittently during routine retrieval of the drill string for 

25 replacement of the drill bit. Sensors responsive to temperature, pressure, gamma radiation, neutrons, acoustic 
energy and electromagnetic energy can be used in the MWD system. The primary disadvantage of MWD over 
wireline measurements lies in the slow rates of data transmission and recording when drilling at fast rates of 
penetration. When MWD data are recorded in real time, telemetry to the surface is typically via the drilling 
fluid column since there is no direct electrical connection between the downhole sensors and the surface 

30 equipment. Current drilling fluid telemetry rates are of the order of bits per second while wireline telemetry 
rates over the logging cable are orders of magnitude greater. In addition, MWD data are created during equally 
spaced time intervals, whether telemetered directly to the surface or recorded downhole for subsequent sur- 
face retrieval and processing. To display MWD data as a f unctton of depth and utilize well established depth 
based signal processing and analysis technology implemented at surface located equipment, MWD data must 

35 be accurately transformed from the dynamic time domain to the depth domain. The time to depth domain trans- 
formation is not a simple function since the penetration rate of the drill bit is neither constant nor predictable. 
In addition, transformation can result in gaps in the depth based data resulting from slow MWD data acquisition 
rates. When drilling at slow rates of penetration, MWD data can be acquired at a much finer axial resolution 
than wireline data. Existing wireline processing techniques are often inappropriate for MWD data. The present 

'M inventive disclosure pertains to new processing techniques which have been developed for the correct con- 
veraton of normally evenly spaced time based data to evenly spaced depth based data which maintain the max- 
imum amount of information from the raw data as well as define the quality of the data. 

Seismology involves the measurement of reflected or refracted acoustic waves in the time domain and 
the conversion or "migration" of these data to the depth domain. Acoustic energy is imparted at the surface of 

45 the earth and reflections or refractions from subterranean features exhibiting contrasts in acoustic impedance 
(the product of acoustic velocity and density) are measured at the earth surface in the time domain. Time do- 
main measurements are then converted to depth domain thereby giving the depth of the perturbing subterra- 
nean features. Depth conversions are only as accurate as the knowledge of acoustic Impedance of all inter- 
vening geological strata. U.S. Patent No. 5,229,940 to Shein S. Wang and David W. Bell teaches a method of 

50 generating three dimensional seismic images of earth formations from a series of conventional two dimen- 
sional seismic traces or lines. The methods are directed solely to seismic technology wherein no well boreholes 
are involved. Seismic processing involves the conversion of reflected or refracted acoustic waves measured 
as a function of time into amplitude as a function of depth using offset velocity measurenrtents and assump- 
tions. The current invention involves the conversion of borehole environs data measured as a function of time 

55 into corresponding data as a f unctk>n of depth using well defined time to depth convereion steps. 

The analysis of returned drilling fluid in a well borehole drilling operation is commonly known as "mud" 
logging. The returned drilling fluid is analyzed for traces of hydrocarbons including gas, and the returned drill 
cuttings are analyzed to determine the porosity, lithology and permeability of the strata being penetrated by 
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the drill bit. Ruid and cuttings are analyzed at the surface of the earth and must t>e assigned a depth of origin 
in order to correlate with conventional wireline logs. MWD logs and seismic data. U.S. Patent No. 5,237,539 
to Thomas H. Selman teaches methods of determining the "lag** time for the returned drilling fluid to travel 
from the drill bit to the surface. This, when correlated with the known depth of the drill bit, allows one to assign 
5 a depth of origin corresponding to each volumetric increment of returned drilling fluid. The subject patent also 
teaches depth correlation of lag time parameters with parameters measured in real time including drill bit depth, 
rate of penetration, weight on bit drill string revolutions per minute and drilling fluid pump pressure. Real time 
and lag time parameters plotted as a function of depth constitute a conventional "mud" log. All of the teachings 
of the subject patent involve measurements at the surface wherein no downhole sensors are employed. All 
10 data are acquired during equal depth intervals rather than during equal time intervals. The patent also teaches 
depth correlation of real and lag time parameters using a simple linear depth shift of the data. There Is no teach- 
ing of quality indicators reflecting the accuracy or reliability of the lag time conversions. 

A primary object of the current Invention is a more correct conversion of MWD data normally acquired at 
evenly spaced time intervals, but unevenly spaced depth Intervals, into evenly spaced depth Intervals so that 
15 readily available methods of data processing, analysis and presentation can be utilized. 

Another objective of the present Invention Is the filling of any gaps In the data resulting from slow data 
acquisition rates when the data are converted to evenly spaced depth intervals. This process involves inter- 
polation of adjacent non-zero data points. 

Another feature of the present invention is the tracking and display of the quality of the data as a function 
20 of depth when the depth sample intervals vary with drill bit penetration rates. Quality parameters include the 
signal to noise ratio, the depth resolution and statistical precision of the data and resulting parameters com- 
puted from the measured data. 

Another aspect of the present Inventbn is the enhancement of raw data resolution when raw data are both 
precise and finely sampled. 

25 Another benefit of the present invention is the proper combination of data from multiple sensors disp>osed 

within the vicinity of the drill bit. Such sensors are usually spaced at different offset axial locations along the 
drill collar and exhibit different measure reference points and axial resolutions. When parameters of interest 
such as formation density, formation porosity, formation resistivity and the like are computed from the com- 
bined responses of two or more sensors, data are shifted to a common measure point and axial resoiutton of 

30 all sensors are matched. 

Still another aspect of the present invention is the alignment of measure points and the matching of axial 
resolutions of two or more parametere measured with a single sensor. Examples are apparent formation den- 
sity and apparent formation photoelectric factor measured with a single sensor using scattered and absorption 
gamma ray techniques. It is well known in the art that these measurements exhibit different effective measure 

36 points and axial resolutions although gamma radlatton is measured with a single sensor. 

An output of the present invention is the tracking of the quality of parameters such as density, porosity, 
resistivity and the like computed from the response of multiple sensors, and the processing applied to each 
response. 

40 Basic Processing 

Transformation of data acquired at evenly or unevenly spaced time Intervals, giving rise to unevenly spaced 
depth intervals, into evenly spaced depth intervals comprises the steps of (a) passing data through a pref liter 
process to place indlvkiual raw data points and groups of points, measured during evenly spaced time intervals, 

45 onto points of a corresponding fine uniformly spaced depth grid, (b) an interpolation process to fill those grid 
points not occupied by data, and (c) then the step of a second filtering and decimatk)n process to put the data 
on an evenly spaced coarse depth grid that is required for normal log analysis. 

The pref liter is a center weighted filter of width 21, where / is the spacing of the fine grid. The center of 
this filter is moved in increments of / along the depth axis containing the unevenly spaced data. If an isolated 

50 data point Is within a distance of tt2 of the center of the filter, the data point is placed at the grid point on which 
the filter is centered. If more than one data point falls within the filter, the data points are weighted by the filter 
shape, points closer to the center of the filter being weighted heavier, and a weighted average data point is 
computed. The weighted average data point is then placed at the filter center depth. 

As a result of the pref iltering process, some grid points contain data while others are empty. Those empty 

55 grid positions are filled by interpolation between adjacent points containing data. As an example, if one point 
is empty and the points immediately above and below contain data, a simple linear Interpolation can be used 
and the empty point is filled with the average of the two adjacent data points, ideally, it would be better to 
avoid interpolation by decreasing the acquisltton time Intervals of the measuring process, thereby decreasing 
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the average depth intervals of the data. Practical problems of processing, data acquisition rates and nuclear 
statistics, if applicable, associated with the measurement lead to lower limits on time acquisition. 

The final step is the decimation of the fine grid data onto a desired coarse grid whose spacing is a multiple 
of the fine grid. Before decimation can occur, a low pass filter is applied to the data on the fine grid to prevent 

5 aliasing. The cutoff frequency, fc, of the filter should be the equal to the Nyqulst frequency (wave number) of 
the coarse grid. That is, fc = 1/2L where L is the spacing of the coarse grid. Adeconvolution process is applied 
to the data to optimize the axial resolution of the measurement. After filtering and decon volution, only those 
points on the fine grid a distance L apart are kept with the remaining points being discarded. As an example, 
if the fine grid has a spacing of 3/4" and the coarse grid has a spacing of 3", only every fourth point of the fine 

10 grid is retained. 

Additional Processing 

For cases where two or more axially spaced sensors are used or more than one parameter is measured 
15 with a single sensor, additional processing is applied to combine the measurements. In both cases, the meas- 
ure points of each measurement may be offset with respect to. each other, in addition, axial resolutions of each 
sensor may differ. As a result, combining measurements requires both depth shifting and resolution matching 
of one measurement with respect to another. 

Depth shifting to a common measurement point can be done on the raw or evenly spaced data. If the raw 
20 data are shifted, the shifted data are then processed as previously described where shifted data are placed 
on the same evenly spaced grid as defined for the reference sensor. If the sensor spacings are an integer mul- 
tiple of the spacing of either the fine or the coarse grid, depth shifting is accomplished with a simple integer 
shift of either the fine or the coarse grid. 

After depth shifting, the data from one or more additional sensors are again filtered to match the depth 
25 resolution of the reference sensor. Resolution matching filters are determined from the Impulse response func- 
tions of the sensors. After depth resolutions are matched, the data from all sensors are then combined to com- 
pute the parameters of interest as a function of depth. 

Quality Control 

30 

As previously mentioned, the quality of the data changes As previously mentioned, the quality of the data 
changes due to variations In drilling rate In depth transformed spacings. The quality of the date Is not generally 
obvious from the log of the processed parameters of Interest. To enhance the value of the log, the quality In- 
dicators (a) precision (repeatebility) (b) depth or axial resolution of the measurements and (c) the signal to noise 

35 are tracked. In addition, this information can be used to control an adaptive filter, with filter characteristics gov- 
erned by some criteria of data quality. Such an adaptive filter may be used to partially deconvolve sensor re- 
sponse, or to reduce noise to some predetermined level. 

The observed error, and statistical error if applicable, of each raw data point can be computed as will be 
deteiled later. This error is tracked when the date are placed onto the fine grid. If several closely spaced points 

40 are averaged together, the error of the averaged point is computed. Also tracked with each point on the fine 
grid is the spacing between raw date points. This spacing represents an "irreducible" resolution length. That 
is, vertical deteils of the penetrated earth formation finer than the irreducible resolution can not be resolved. 

When interpolated points are generated, appropriate spacing numbers and precision numbers are as- 
signed. After the anti-aliasing or other filter or deconvolution is applied to the date, the precision errors and 

45 resolution lengths of the data are recomputed and tracked. Both precision errors and resolution lengths are 
presented on the log along with the petrophyslcal parameters of interest thereby providing the log analyst with 
valuable indicators of the quality of the measured data. 

So that the manner in which the above recited features, advantages and objectives of the present invention 
are attained and can be understood in detail, a more particular description of the inventton, briefly summarized 

50 above, may be had by reference to the embodiments thereof which are illustrated in the appended drawings. 
It is to be noted, however, that the appended drawings illustrate only typical embodiments of this invention 
and are therefore not to be considered limiting of its scope, for the invention may admit to other equally effective 
embodiments. 

Fig. 1 1s a generalized lllustratton of a drilling system employing MWD for determination of earth formation 
55 and borehole parameters; 

Fig. 2 depicts a flow chart for data processing as detailed In the current disclosure; 

Fig. 3 illustrates a comparison of apparent formation density as computed from the response of a short 

axially spaced sensor and a long axially spaced detector using the raw data as recorded in equal time in- 
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tervals and using data transformed to equal depth intervals using the methods of the invention; 
Fig. 4 illustrates a log of apparent formation density as computed from the response of a short axially 
spaced sensor and a long axially spaced sensor measure point alignment and sensor resolution matching 
using methods of the invention; 
6 Fig. 5 illustrates a log of formation density which has been compensated for adverse effects of the borehole 

by combining the apparent formation densities computed from the responses of the short spaced and long 
spaced sensors; and 

Fig. 6 illustrates logs of the quality of the apparent density measurements made by the short axially spaced 
and long axially spaced sensors. 

10 Data collected from downhole sensors in MWD operations are normally collected during evenly spaced 

time intervals. These data are either telemetered to the surface in real time for processing or are recorded 
downhole for retrieval and processing when the drill string is pulled from the well bore to change the drill bit. 
Data can also be collected during unevenly spaced time Intervals. This requires the data to be normalized to 
a unit time interval before further processing as detailed betow. Data acquisition may also be intermittent An 

15 example would be termination of data collection upon interim cessation of drilling. Again, data must be nor- 
malized to a unit time interval before further processing. The preferred embodiment is the collection of data 
in evenly spaced time intervals. 

Analysts are interested in formation parametric measurements as a function of depth within the well bore. 
This permits the analyst to correlate MWD logs with conventional wireline logs of offset wells, seismic data, 

20 geological maps and other information used to evaluate the well being drilled. This information regarding ad- 
jacent wells will always have depth related characteristics. In addition, data processing and interpretation tech- 
niques which have been developed over the yeare are, in general, depth based resulting from their origins in 
wireline logging and seismic processing, it is advantageous, therefore, to convert MWD data to the depth do- 
main rather than time base. The conversion of MWD data from the time to depth domain is not a simple process. 

25 This is due to the fact that the penetration rate of the drill bit is not constant in time. Time to depth domain 
conversion of MWD data, in addition to treating variable bit penetration rate, must also preserve the precision 
and axial or vertical resolution of the various downhole measurements. In measurements involving nuclear 
processes, optimization of the statistical precision of the measurements Is especially critical. Finally, if the re- 
sponses of two or more sensors are used to compute the parameters of Interest, vertical resolution of the sen- 

30 sors must be matched to avoid spurious signals while still optimizing resolution and signal to noise ratio. The 
current invention is able to meet all of these key criteria. 

Fig. 1 illustrated a typical drilling system employing MWD apparatus. The MWD tool 10 is mounted on the 
drill string 12 in the immediatevicinlty of the drill bit 13. A standard drilling rig (not shown) is employed for control 
of the drill string In the borehole 14 within an earth formation 18. The invention as disclosed Is applicable to 

36 data from one or more sensore within the MWD tool 1 0. For purposes of Illustrating the preferred embodiment, 
a short axially spaced sensor 15, a second long axially spaced sensor 16, and a radiation source 19 are de- 
picted. The sensors 15 and 16 are powered with suitable power sources (not shown) within the tool 10. Data 
are either stored in a recorder (not shown) within the tool 1 0 for later retrieval or are telemetered by a transmitter 
17 to the surface to a receiver 11 for real time processing by a central processing unit (CPU) 22. After proc- 

40 essing using methods to be detailed, the output parameters of interest are passed to a recorder 24 wherein 
the log 25 of the parameters of Interest as a function of depth is produced. 

Aflow chart of the data processing algorithms associated with the current Invention Is shown In Fig. 2. In 
the preferred embodiment, raw data f| from a sensor which is acquired at evenly spaced time intervals, but 
unevenly spaced depth intervals. Corresponding depth intervals are assigned to the raw data fi at step 29. 

45 These data are next placed onto a corresponding fine, uniformly spaced depth grid with spacing / between 
points. Note that being unevenly spaced in depth, the data set is not aligned with the evenly spaced grid points 
at this stage of the processing. This step is illustrated at 30 of the flow chart. The error e| associated with each 
corresponding value of f| is next computed in step 32. For nuclear data, the well known Poisson statistical error 
associated with the counts is computed at 32. The spacing ri between each point f| on the fine grid Is computed 

50 at the step 34. Both e| and n are tracked throughout processing and eventually result In quality logs which will 
give the analyst an Indication of the error and vertical resolution of the quantities of Interest. 

The data are pref lltered at the step 34 to place Individual raw data points and groups of raw data points, 
measured during equally spaced time intervals, onto the points of the fine equally spaced depth grid. The pre- 
ferred pref liter is a center weighted filter of length 2/. The filter is moved in increments of 1 along the fine grid. 

55 If an isolated data point is within a distance of 1/2 of the center of the filter, the data point is placed at the depth 
of the center of the filter. >f more than one data point falls within the defined weighted filter, the data points 
are weighted by the filter shape, points cfoser to the center of the filter being weighted heavier, and the weighed 
average is computed. Stated mathematteally: 
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j 

where f^j is the value of the pref iltered data point assigned to ghd point i of the fine grid, w^j are the weighting 
functions and J denotes the raw data points falling within the filter centered at grid point i. If one assumes that 
all values of AU| the spacings between raw data points, are equal, then equation (1) reduces to 

(2) 

J 



The intervals AU| are not equal as previously discussed. In practice, however. If the length of the filter is small 
when compared to the actual measured variation of the data, then the assumption Introduces negligible error. 
The error ea,i and resolution r^ , associated with pref iltered data fa,i are computed at the step 36. 

As a result of pref iltering, some points of the fine grid contain no data. Each grid point is tested at the step 
38. Those empty grid positions are filled by mathematical interpolation between the nearest adjacent non-zero 
grid points fa.(.j and fa.|.t.k step 40. One such Interpolation scheme is given In equation (3); 

^ ^ - /)] + mk - i)i 

where i denotes the fine grid point being filled by interpolation, and the Indices l-J and i+k denote the nearest 
adjacent grid points containing data on either side of grid point i. Associated error and resolution are computed 
at the step 42. An observed or "noise" error Is also estimated from the variance of a line fit to all non-interpolated 
points within the detector response length. This gives the noise level, with the rationale being that negligible 
changes due to the formation will occur within the response length of the detector. Poisson statistical error is 
also computed for interpolated and non-interpolated points. Poisson statistical error for interpolated data points 
are obtained by well known statistksai error propagation techniques. The systematic or noise error is used for 
Oaj unless the noise ievei is less than the statistical error, in this situation, eaj is set to equal the statistlcai 
error. Many techniques can be used to calculate the resolution parameter ra,i ranging from a simple traddng 
of spacing between raw data points to more complex expresstons involving sensor response functions, asso- 
ciated error and sample rate. 

At this juncture, the fine grid of spacing / has been filled with data. Thefine grid is then decimated onto a 
more conventional coarse grid once the fine grid data have been processed with anti-aliasing filter, see the 
step 44. Signal processing theory states that the Nyquist frequency fc is the maximum frequency supportable 
on a grid of spacing L and further that fg = 1/2L An anti-aliasing filter that exhibits a sharp cut-off at the Nyquist 
frequency Is applied to further filter the data and obtain fb,i. The function 8in(xrr)/x which Is centered on a / ^ 
3/4" grid of in the range of 20 to 60 points in the preferred embodiment The number of points can be varied 
without significant loss of information so long as the filter length is longer than the correlation length of the 
data. The Fourier transform of a function of the form sin(ax)/x yields a step function of amplitude n when cd is 
less than a and greater than -a. The length constant T can therefore be set to give the desired cut-off frequency 
CO = 2 fcfc where in the preferred embodiment, fc is the Nyquist frequency of a coarse grid with spacing L= 3.0". 
Therefore T = 1/w = V{2icf^. Since fc = 1/2L, T=\Jn = 0.995" for the preferred embodiment A deconvdution 
process is also applied to the data at step 44 in order to enhance the axial resolutton of the measurement 

After anti-aliasing filtering and deconvoiution, the associated errors, eu, and resolutions, r^^ are computed 
at the step 46. 

Continuing with the description, the step 48 is the decimation of the fine grid / onto the coarse grid L where 
L is a multiple of/. After filtering, only those points on thefine grid with a spacing of Lare kept with the remaining 
points being discarded. In the preferred embodiment where / = 3/4" and L = 3", only every fourth point on the 
fine grkl is kept. This decimated data set is now designated as with associated error and resolution values 
of Em and respectively. 

Going further with the description, the step 50 queries the presence of more than one sensor or if more 
than one parameter is measured using the same sensor, if additional sensors or multiple measurements from 
a single sensor are used, previously explained steps 30 through 48 are repeated for the second or subsequent 
set of raw data. Multiple sensors axiaUy spaced within tool 1 0 or multiple parameter measurements from a sin- 
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gle sensor will usually exhibit different effective measure points. When Gon>blnlng responses of two or more 
data sets to obtain the parameter of Interest, it is necessary to shift the second and subsequent data sets so 
that the effective measure point(s) are aligned with the measure point of a designated reference sensor or 
measurement. Depth shifting can be done on raw data or on evenly spaced data resulting from filtering at step 

5 34 or step 44. If raw data from a subsequent sensor are shifted, processing steps 30 through 48 are applied 
to the shifted data set. If the effective sensor spadngs are Integer multiples of either the fine grid spacing or 
the coarse grid spacing, depth shifting is accomplished with a simple integer shift of either the fine or the 
coarse grid. The latter event Is depicted as step 54 in Fig. 2. 

One additional process must be applied to data from multiple sensors or multiple measurements from a 

10 single sensor before they are combined to compute parameters of interest. The vertical responses or vertical 
resolution of multiple sensors are not necessarily identical and, in the majority of applications, are significantly 
different. Combining multiple responses across any geological boundary penetrated by the drill bit will give 
unrealistic results unless the axial resolutions are matched. The procedure is to match the higher resolution 
measurement to the lower resolution measurement. In the case of applications involving Induced nuclear meas- 

15 urements, the higher resolution sensor closer to the source of radiation is matched to the sensor spaced ferther 
from the source of radiation. Resolution matching is accomplished at step 56 by applying a convolution filter 
to the high resolution sensor. The convolution filter is defined as a function, when convolved with the response 
function of the high resolution sensor, that yields the response of the low resolution sensor. In frequency space, 
this Is a convolution equated to a product. The frequency components of the unknown resolution matching f unc- 

20 tlon are calculated by dividing the frequency components of the low resolution sensor impulse response by 
the frequency components of the high resolution impulse response. 

At this point In the processing algorithm, the data which were originally measured In equally spaced time 
Intervals have been converted to evenly spaced depth Intervals. In the case of multiple sensors, data have 
been depth shifted and the resolution(s) of the sensors have been matched. The parameter of interest, des- 

25 ignated as Pm, is computed at the step 60 from the response of one sensor or by combining the response of 
two or more sensors. Error associated with each data point has been computed and propagated through the 
various processing steps so that error Em Is assigned to each corresponding value of Pm. Spacing between 
raw data points has also been tracked through the various processing steps allowing a vertical resolution para- 
meter Rjn to be assigned to each corresponding value of P^. Pm and the quality parameters E^n and Rgn are 

30 recorded as a function of depth by the recorder 24 at step 62. 

As stated previously, tracking of the quality indicator error and resolution is a valuable aid to the analyst 
in interpreting MWD logs. Quality indicators are also valuable means for driving adaptive filters. Error is not 
only a function of sensor design and sample rate but it is also a function of the physical characteristics of the 
formation. Using induced nuclear measurements as an example, measured counting rates control the statlst- 

35 leal precision of the measurement and are also a function of the characteristics of the formation. In addition, 
it Is well known that the response of a sensor moving across a bed boundary In an Induced nudear measure- 
ment is altered by the contrast of the adjacent beds. Stated in another way, the contrast at bed boundaries 
can impact vertical resolution. In addition, formations can vary the imputed measure points of sensors in in- 
duced nudear applications. Quality indicators as disdosed are, therefore, ideally suited to adaptively vary filter 

40 characteristics depending upon actual downhole conditions encountered. Adaptive filtering is used to maxi- 
mize statistical precision and axial resolution, and is especially applicable in the preferred emt)odiment at the 
depth shift step 54 and resolution nnatchlng step 56 of the processing algorithm. 

Applications 

45 

A formation density measurement, compensated for the adverse effects of the borehole, will be used to 
further illustrate the preferred embodiment of the invention. The measurement equipment involves two gamma 
ray sensors or detectors axially spaced from a source of gamma radiation mounted in the immediate vicinity 
of the detectors. These detectors are known as the short spaced and long spaced detectors to those skilled 

50 in the art. The physical basis of the gamnra-gamma dual detector compensated density measurement was first 
presented in the paper *The Dual Spacing Formatton Density Log", J. S. Wahl, J. Tittman, and C. W. Johnstone, 
Journal of Petroleum Technology, December, 1964 and is hereby referenced, it is again emphasized that the 
current invention is not limited to the dual detector compensated density measurement or to nuclear measure- 
ments in general, but it is applicable to measurements made by one or more sensors of any type in the time 

55 domain. 

Fig. 3a illustrates the plot 63 of apparent formation density computed from raw data from the short spaced 
detector plotted as a function of depth In the well borehole in meters. These data are compared with the same 
raw data processed through the declmatton step 48 and computattonal step 60. Fig. 3b Illustrates a similar 
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plot from the long spaced detector. It is apparent that both the processed short and long spaced density meas- 
urements closely track the raw data. It should be noted, however, that the short spaced detector exhibits su- 
perior vertical resolution as Is exhibited by the relatively small full- width- half-maxima 65 compared with 66 of 
the major density excursions. It should also be noted that the measure point of the short spaced detector is 

5 "deeper" than that of the long spaced detector. This Is due to the fact that in the example the short spaced 
detector is closer to the source of gamma radiation which, in turn, is doser to the drill bit 13. If one attempted 
to combine the data depicted in Figs. 3a and 3b to obtain a density computation.or to use other combinations 
of multiple parameters with depth resolution mismatches, serious error would be introduced as a result of re- 
solution mismatch and data shift This error would be especially pronounced in zones where density is rapidly 

10 changing. 

Fig. 4 shows apparent fomnation density 70 and 72 computed from the response of the long and short 
spaced detectors, respectively, as a function of depth after the depth shifting step 54 and resolution matching 
step 56. It is apparent that the two sensors are on depth and that resolution has been matched as illustrated 
by the near overlays of density excursions at 45B5.3 m, 45B6.0 m and 45B6.7 m. 

15 Fig. 5 shows a borehole compensated formation density computed from the apparent formation density 

as measured with the short and long spaced detectors using computational concepts previously referenced. 
It is apparent that the compensated density curve exhibits no unusual excursions such as "spikes" or "horns" 
at bed boundaries indicating proper depth shifting and resolutk>n matching. 

Fig. 6 shows a quality logs of corresponding error Indicators and resolution indicators R^, again plotted 

20 as a function of depth within the well borehole in meters. Fig. 6a depicts apparent formation density 83 as meas- 
ured by the short spaced detector, while the curves 81 and 85 define the one observed standard deviation 
uncertainty limits associated with the parameter of Interest. Using previously defined nomenclature, the curve 
81 represents +Ejn and the curve 85 represents - E,n. Fig. 5b is the corresponding plot of associated 
error Em curve 89 and resolution 87 associated with the apparent density measurement 83 from the short 

26 spaced detector. Fig. 6c illustrates the one standard observed uncertainty limits 91 and 95 of the apparent 
formation bulk density 93 measured by the long spaced detector, while Fig. 6d illustrates the corresponding 
error 99 and resolution 97 associated with the density measurement 93. It is again emphasized that error and 
resolution calculations can also be made for the compensated bulk density as well as apparent densities cal- 
culated from the long and short spaced detectors. 

30 Those skilled in the art will recognize modifications and refinements to the preferred embodiment, includ- 

ing adaptive filtering, variations of fixed filter types, and variations in the number, types and positk>ns of the 
sensors. Such enhancements and modifications are within the scope of the present invention as identified in 
the following claims. 

36 

Claims 

1. A method for translating raw data acquired from a MWD system operating in a well borehole wherein the 
raw data are acquired at spaced time intervals and are translated onto corresponding evenly spaced depth 

40 Intervals, the method comprising the steps of; 

(a) acquiring raw data in spaced time intervals using a sensor in a well boretioie and comprising a part 
of a MWD system; 

(b) assigning to said raw data, depth values corresponding to the locations of acquisition; 

(c) placing said raw data onto a first evenly sp3ced depth grid; 

45 (d) positioning data onto points of said first depth grid by filtering said raw data using a first depth filter 

thereby creating a second data set; 

(e) fbr selected points of said first evenly spaced depth grid not containing data, providing data for saM 
points by interpolating adjacent points of said first evenly spaced depth grid containing data; 

(f) applying a second filter, a deconvolutton process or both to said second data set creating a third 
so data set; 

(g) placing said third data set onto a second evenly spaced depth grid with spacing greater than said 
first evenly spaced depth grid; and 

(h) retaining elements of said third data set which fall upon the grid points of said second evenly spaced 
depth grid thereby creating a fourth data set. 

55 

2. The method of claim 1 wherein the errors associated with said raw data acquired at spaced time Intervals 
are tracked through each step in the translatk>n of said raw data to said fburth data set 
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3. The method of claim 2 wherein the axial resolution associated with said raw data collected at spaced ti me 
Intervals are tracked through each step of the translation of said raw data to said fourth data set. 

4. The method of claim 2 wherein said traclced errors are used to adaptively adjust the first filter, the second 
filter, the deconvolution process or any combination thereof. 

5. The method of daim 3 wherein said tracked resolutions adaptively adjust the first filter, the second filter, 
the deconvolution process or any combination thereof. 

6. The method of daim 3 wherein said tracked errors and said tracked resolutions are used to adaptively 
adjust the first filter, the second filter, the deconvolution process or any comblnatton thereof. 

7. A method for translating two or more raw data sets acquired with a MWD system operating In a well bore- 
hole where the rawdata are acquired atspaced time intervals and are translated onto corresponding even- 
ly spaced depth intervals, the method comprising the steps of; 

(a) acquiring, in timed intervals, sets of raw data from one or more sensors mounted within said MWD 
system; 

(b) assigning to saki raw data In each set depth values corresponding to the locations of acquisition; 

(c) placing saki raw data from each said set onto a first evenly spaced depth grki; 

(d) positioning data onto f)olnts of said first evenly spaced depth grid by filtering each said raw data 
set using a center weighted first filter thereby creating second data sets; 

(e) for any point of said first evenly spaced depth grid not containing data for a given second data set, 
filling said point for each given said second data set by interpolating adjacent points containing data 
within said given second set; 

(f) applying a second center weighted filter, a deconvolution process or both to each said second data 
set thereby removing any aliasing and thereby creating third data sets; 

(g) placing each of said third data sets onto a second evenly spaced depth grkl with spacing an integer 
multiple of said first evenly spaced depth grkI; and 

(h) retaining elements of each of said third data sets which fall upon the grid points of said second evenly 
spaced depth grid thereby creating fourth data sets. 

8. The method of claim 7 wherein the errors associated with each said raw data set are tracked through each 
step in the translation of raw data sets to said fourth data sets. 

9. The method of daim 8 wherein the axial resolutions associated with each said raw data set are tracked 
through each step of the translation of said raw data sets to said fourth data sets. 

1 0. The method of claim 8 wherein said tracked errors are used to adaptively adjust the first filter, the second 
filter, the deconvolution process or any combinatk>n thereof. 

11. The method of daim 9 wherein saM tracked axial resolutk>ns are used to adaptively adjust the f Iret filter, 
the second filter, the deconvolution process or any combination thereof. 

12. The method of daim 9 wherein said tracked errors and said tracked axial resolutions are used to adap- 
tively adjust the first filter, the second filter, the deconvolution process or any combinatton thereof. 

13. The method of daim 7 wherein said raw data sets exhibit different axial measure points and axial reso- 
lutions and: 

(a) said fourth data sets are shifted to a common axial measure point thereby forming fifth data sets; 

(b) saki fifth data sets are processed to match axial resolutions of the sensore thereby forming sixth 
data sets; 

(c) said sixth data sets are combined to yield computed parameters of Interests from the borehole en- 
virons and earth formatk>ns traversed by said MWD system; 

(d) errors and axial resoiutk>ns associated with each data set are tracked from said first raw data sets 
through said sixth data sets and are used to determine parametric errore and parametric axial resolu- 
tions of said computed parameters of interest; and 

(e) said parameters of interest and said parametric errore and parametric axial resolutions are recorded 
as a function of depth within sakI well borehole. 
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14. The method of claim 7 wherein the raw data are acquired in evenly spaced time intervals. 

15. A system for acquiring and translating one or more raw data sets from at least one sensor In a MWD system 
operating In a borehole wherein the raw data are acquired at spaced time Intervals and are translated 
onto corresponding evenly spaced depth Intervals, comprising: 

(a) sensor means for a MWD system for acquiring said raw data from a well borehole; 

(b) means for receiving said raw data from said sensor or sensors; 

(c) means of assigning a depth value to each measured raw data value corresponding to the location 
of acquisition; 

(d) means for positioning each of said raw data sets onto a first evenly spaced depth grid; 

(e) means for positioning data onto points of said first evenly spaced depth grid by filtering said raw 
data sets using a center weighted first filter thereby creating second data sets; 

(0 for any point of said first evenly spaced depth grid not containing data for a given said second data 
set, means of filling said point for each given second data set by interpolating adjacent points containing 
data within said second data set; 

(g) means for applying a deconvolutlon process, a second center weighted filter, or both to each said 
second data set thereby removing any aliasing and thereby creating third data sets; 

(h) means for placing said third data sets onto a second evenly spaced depth grid with the spacing of 
said evenly spaced depth grid being an integer multiple of said first evenly spaced depth grid; 

(i) means for retaining elements of each said third data sets which fall upon a grid point of said second 
evenly spaced depth grid thereby creating a fourth data sets; and 

(j) means for recording said fourth data sets. 

16. The system of claim 1 5 further comprising means for tracking errors associated with said raw data sets 
through each step of the translation of said raw data sets to said fourth data sets. 

17. The system of claim 1 6 further comprising means for tracking axial resolutions associated with said raw 
data sets through each step of the translation of said raw data to said fourth data set. 

18. The system of claim 16 further comprising means for using saki tracked errors to adaptively adjust the 
first filter, the second filter, the deconvolutbn process or any combination thereof for each data set. 

19. The system of dalm 17 further comprising means for using said tracked axial resolutbns to adaptively 
adjust the first filter, the second filter, the deconvofution process or any combinatton thereof for each data 
set. 

20. The system of claim 1 7 further comprising means for using sakJ tracked errors and said tracked axial re- 
solutions to adaptively adjust the first filter, the second filter, the deconvolutlon process or any combin- 
ation thereof for each data set 

21. The system of claim 1 7 wherein two or more sensors are spaced at different axial focations within said 
MWD system and further comprising: 

(a) means of adjusting sakJ first filter, said second filter, said deconvolutlon process or any combination 
thereof based upon the values of sakJ tracked errors and said tracked resolutions for each said sensor; 

(b) means for shifting said fourth data sets to a common axial measure point thereby forming fifth data 
sets corresponding to the response of each said sensor at common measure points; 

(c) means for processing said fifth data sets to match axial resolutions of said sensore thereby forming 
a sixth data set corresponding to the response of each said sensor; 

(d) means for combining said sixth data sets from each said sensor to yield computed parameters of 
interest of earth formations and borehole environs travereed by the MWD system; 

(e) means for tracking sakI errors and said axial resolutions associated with said raw data sets to said 
sixth data sets; 

(f) means for converting said errors and said axial resolutions associated with said sixth data sets into 
parametric errore and parametric resolutions of said computed earth formation and borehole environs 
parameters of interest; and 

(g) means for recording said parameters of interest and said parametric errors and sakI parametric axial 
resolutk>ns as a functk>n of depth of sakI common measure point of the MWD system within said well 
borehole. 
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22. The system of claim 1 5 wherein raw data are acquired at evenly spaced time Intervals. 

23. The system of claim 15 wherein said sensors acquire raw data induced by a source of radiation within 
the MWD system. 

24. The system of claim 15 wherein said sensors acquire raw data induced by a source of acoustic energy 
within the MWD system. 

25. The system of claim 1 5 wherein said sensors acquire raw data induced by a source of electromagnetic 
data within the MWD system. 

26. The system of claim 15 wherein said sensors acquire raw data resulting from naturally occurring radio- 
active material within earth formations traversed by the MWD system. 

27. A method of translating raw data acquired from a MWD system operating in a well borehole wherein the 
raw data are acquired at spaced time intervals and are translated onto corresponding evenly spaced depth 
intervals, the method comprising the steps of: 

(a) acquiring raw data in spaced time intervals using a sensor in a well borehole and comprising a part 
of a MWD system: 

(b) assigning to said raw data, depth values corresponding to the locations of acquisition; 

(c) placing said raw data onto an evenly spaced depth grid; 

(d) positioning data onto points of said depth grid by filtering said raw data using a firet depth filter there- 
by creating a second data set; 

(e) for selected points of said evenly spaced depth grid not containing data, providing data for said 
points by interpolating adjacent points of said evenly spaced depth grid containing data; and 

(f) applying a second filter to said second data set thereby creating a third data set. 

28. The method of claim 27 wherein the errors associated with said raw data acqu ired at spaced time intervals 
are tracked through each step in the translation of said raw data to said second data set 

29. The method of claim 28 wherein the axial resolution associated with said raw data collected at spaced 
time intervals are tracked through each step of the translation of said raw data to sakJ second data set 

30. The method of claim 28 wherein said tracked errors are used to adaptively adjust the first, second or both 
of said depth filters. 

31. The method of daim 29 wherein said tracked resolutions are used to adaptively adjust the first, second 
or both of said depth filters. 

32. The method of claim 29 wherein said tracked errors and said tracked resolutions are used to adaptively 
adjust the first, second or both of said depth filtere. 

33. A system for acquiring and translating raw data from at least one sensor in a MWD system operating in 
a well borehole wherein the raw data are acquired at spaced time intervals and are translated onto cor- 
responding spaced depth intervals, comprising: 

(a) sensor means for a MWD system for acquiring said raw data from a well borehole; 

(b) means for receiving raw data from said sensor or sensors; 

(c) means for assigning a depth value to each measured raw data value corresponding to location of 
acquisitk>n; 

(d) means of positioning raw data onto an evenly spaced grid; 

(e) means of positioning data onto points of said evenly spaced grid by filtering said raw data using a 
first center weighted depth filter thereby creating a second data set; 

(f) for any point of said evenly spaced depth grid not containing data for said second data set, means 
of filling said point by interpolating adjacent points containing data within said second set; 

(g) means for applying a second center weighted depth filter to said second data set thereby removing 
any aliasing and thereby creating a third data set; and 

(h) means of recording said third data set 

34. The system of claim 33 further comprising means of tracking errora and resolutions associated with said 
raw data sets through each step of the translation of said raw data to sakJ third data set 
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35. The system of claim 34 further comprising means of using said tracked errors, tracked resolutions or both 
to adaptively adjust the first second or both of said depth filters. 
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